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Abstract. Several experiments in the context of ladder materials have recently shown that the study of 
simple models of anisotropic ladders (i.e. with different couplings along legs and rungs) is important for the 
understanding of these compounds. In this paper Exact Diagonalization studies of the one-band Hubbard 
and t — J models are reported for a variety of densities, couplings, and anisotropy ratios. The emphasis 
is given to the one-particle spectral function A(q, u)) which presents a flat quasiparticle dispersion at the 
chemical potential in some region of parameter space. This is correlated with the existence of strong pairing 
fluctuations, which themselves are correlated with an enhancement of the bulk-extrapolated value for the 
two-hole binding energy as well as with the strength of the spin-gap in the hole-doped system. Part of the 
results for the spectral function are explained using a simple analytical picture valid when the hopping 
along the legs is small. In particular, this picture predicts an insulating state at quarter filling in agreement 
with the metal-insulator transition observed at this special filling for increasing rung couplings. The results 
are compared against previous literature, and in addition pair-pair correlations using extended operators 
are also here reported. 

PACS. 71.27.+a Strongly correlated electron systems - 74.72.-h High-Tc cuprates - 75.40.Mg Numerical 
simulation studies - 79.60.-i Photoemission and photoelectron spectra 

1 Introduction understand the mechanism of superconductivity in strongly 

correlated low-dimensional systems. 

Ladder structures can also be found in other oxides 
The discovery of superconductivity in quasi-two dimen- such as vanadates |. Magnetic susceptibility measure- 
sional (2D) copper-oxide materials has led to a renewed ments 011 M S V 20s have been interpreted in terms of weakly 
interest in the physics of doped Mott-Hubbard insulators cou P led H «senberg ladders. In addition, recent X-rays 



H i and in the interplay between magnetism and supercon- scattering experiments ^ have suggested that NaV 2 
■ ductivity. Recently another class of copper oxide materi- could be considered as a quarter-filled ladder system, 
als based on weakly coupled one dimensional (ID) lad- While the stoichiometric parent compounds of the su- 
ders has been synthetised. Their structure is closely perconducting 2D cuprates are antiferromagnetic Mott in- 
related to the one of the 2D perovskites, namely it con- sulators, the parent insulating ladders exhibit spin liq- 
tains S=l/2 copper spins which are antiferromagnetically uid properties. The existence of a spin gap in a spin- 
coupled along the ladder direction (legs) and along the ladder structure has been first proposed theoretically §j 
rungs through Cu-O-Cu bonds. Recent experimental re- and found experimentally in several even-leg ladder copper- 
suits reporting superconductivity § in the hole-doped oxide systems (such as SrCu 2 3 §J| and LaCu0 2 . 5 0)- 
ladder cuprate Sri4_ :c Ca 2: Cu2404i have clearly established It has been suggested that the spin gap, if robust un- 
that the existence of a superconducting state is not re- der doping, could be responsible for an attractive inter- 
stricted to 2D doped antiferromagnets but it actually ex- action between holes on the same rung §,§. Although 
tends to a wider class of strongly correlated copper-oxide recent experiments [^) suggest that the spin gap disap- 
materials. Thus, the novel ladder compounds offer new pears in hole-doped ladders at the high pressure needed 
perspectives, both for experimentalists and theorists, to to achieve superconductivity, part of the spin excitations 
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are still suppressed as the temperature is lowered in the 
normal state as predicted theoretically . Such a behav- 
ior bears similarities with the pseudogap behavior of the 
underdoped 2D superconducting cuprates. 

It is expected that the formation of hole pairs on the 
rungs can lead to competing superconducting pairing or 
4kp Charge Density Wave (CDW) correlations as e.g. found 
in the weak coupling limit Jljj . For isotropic t — J lad- 
ders (i.e. with equal couplings along legs and rungs), it 
has been established that the d- wave-like superconducting 
pairing correlations [12] are dominant in a large region of 
the phase diagram [Of. Such a state, which is referred to 
as CI SO in the language of the RG analysis (wj, is charac- 
terized by a single gapless charge mode and a gap in the 
spin excitation spectrum |Io| and belongs to the same uni- 
versality class as the Luther-Emery phase of the ID chain. 
Because of the one-dimensional character of the system, 
the superconducting correlations at large distances still 
behave following power laws Jll|jr^ ] . Using finite size scal- 
ing analysis and conformal invariance relations |]l3| , the 
corresponding critical exponents have been computed in 
the case of equal rung and leg couplings (isotropic case). 
Note that a small Josephson tunneling between the lad- 
ders is expected to give a finite superconducting critical 
temperature [|l5) . At finite doping, the spin gap is expected 
to vanish below a small critical J/ 1 ratio jl3|,[l6| . Presum- 
ably, such a transition is associated to an instability of 
the hole pairs gas towards a liquid made out of individual 
holes with two spin and two charge collective modes (C2S2 
phase) [[l6). Whether the disappearance of the spin gap 
observed by NMR experiments in the doped Sri4Cu2404i 
superconducting ladder material under pressure j^] is con- 
nected to such a transition is under much debate. 

Although the isotropic case is the most widely ana- 
lyzed situation in the context of theories for ladders, it has 
now become clear that most of the actual ladder materi- 
als have in fact different leg (J||) and rung (JjJ magnetic 
couplings and/or hopping matrix elements (legs) and 
t± (rungs). Recent neutron scattering experiments |17| ] on 
the insulating ladder Sri4Cu2404i actually suggest a ra- 
tio of — 0.5. A similar anisotropy was in fact also 
predicted previously in the context of the SrCu203 mate- 
rial fis|j . On the other hand, the vanadate ladder NaV205 
apparently corresponds to the opposite limit of strong 
rung couplings with a ratio tj_/i|| — 2 [[l9) which could 
justify a description at quarter-filling in terms of an effec- 
tive ID Heisenberg model 119|J2C 



1. 



A e (q,w) = Im( 



The purpose of the present paper is to investigate spec- 
tral properties and superconducting correlations of anisotropigied electronic states, 
Hubbard and t — J ladders by exact diagonalization meth- 
ods. Finite size scaling analysis is used to obtain physical 
quantities such as the pair binding energy or the spin gap. 
Dynamical correlations functions are also computed using 
a continued fraction method. The focus of the paper will 
be on the role of the ladder anisotropy (regulated by the 
ratios of the rung couplings to the leg couplings) as well 
as on the influence of doping. The anisotropic Hubbard 



ladder is defined as, 

H = t \\Yl ( C i a ;* c *+l,a;a + h.C.) 

where the index a stands for the chain index (— 1,2). 
Anisotropy ratios r a = t±/t^ in the range 0.5 < r a < 
2.5 will be considered. Most of the calculations reported 
below have been carried out in the strong coupling regime 
U/t\\ — 8. Motivated by the doped cuprate and vanadate 
ladders, the studies below are performed in the electron 
density range 0.5 < n < 1. For U/t\\ > 1 and U/t± > 1, 
the low energy spin and charge degrees of freedom can be 
described by the effective anisotropic t — J ladder with 
doubly occupied sites projected out, 

H = J|| ^^(S^Q. • Si + i i0 , — — n^ Q rii + i iQ ,) 

+ j± y^(Sj,i • s i:2 - ^,1^,2) (2) 

i 

+ H (4, a -A+l, a -a + h.C.) 



where c\ a . a = Ci ta ^- a {\ — rii^a) ar e hole Guzwiller pro- 
jected creation operators. The large-t/ limit of the anisotropic 
Hubbard ladder leads to antiferromagnetic exchange cou- 
plings of the form Jp = ^tj 3 /U in the two directions (3 =\\ 
,_L. Hence, for simplicity, the relation J±/J\\ = (tj_/tii) 2 
will be here assumed even outside of the range J±/t± <C 1 
and J\\/tu < 1 of rigorous validity of the equivalence be- 
tween the two models. In the rest of the paper, energies 
will be measured in units of t|| unless specified otherwise. 



2 Single particle spectral function 
2.1 Motivations 

Let us examine first the one-particle spectral function, 

A(q, w) = A e (q, u) + A h (q, u) , (3) 
where A e (q,w) corresponds to the density of the unoccu- 



j q;cr 



ie-H 



>0 ■ ( 4 ) 



and A h (q,u)) corresponds to the density of the occupied 
electronic states (ie unoccupied hole states), 
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Here (•■•) stands for the expectation value in the ground 
state wave function of energy E . The transverse compo- 
nent of the momentum q takes only the two values q^ = 
0,7r corresponding to symmetric or anti-symmetric states 
with respect to the reflection exchanging the two chains. 
Ah(q, oj) and A e (q, u) are of crutial importance since they 
can be directly measured in angular-resolved photoemis- 
sion (ARPES) and inverse photoemission (IPES) spec- 
troscopy experiments, respectively. 

Thus far, the role of the anisotropy ratio t±/tu in dy- 
namical properties of ladders has not been studied in de- 
tail, except at half- filling n = 1. In this case, the spin gap 
is remarkably robust and persists down to arbitrary small 
interchain magnetic coupling J± The single particle 
(and two particles) spectral functions of the Hubbard lad- 
der have been obtained at n = 1 using quantum Monte 
Carlo (QMC) simulations Uj. Working at U = 8, two 
regimes were identified |2^] depending on the magnitude of 
the ratio t±/tu. For instance, increasing t± the half-filled 
Hubbard ladder evolves from a four-band (at small t±/tii) 
to a two-band (at large tj_/tn) insulator. The latter regime 
can be understood from the non-interacting U ~ picture: 
in this case, two heavily weighted bonding (q± = 0) and 
anti-bonding (q± = 7r) bands are separated by an energy 
~ 2t± and the Fermi level lies in between. On the other 
hand, a small fraction of the total spectral weight is found 
in the inverse photoemission spectrum (u> > /i) for q± = 
and in the photoemission spectrum (to < /j) for q± = tt. 
In fact, in the t±/t^ > 1 limit, the spin-spin correlation 
length is very short [^2| and a description in terms of a 
rung Hamiltonian (reviewed in the next section) is appro- 
priate (£|| can then be treated as a small perturbation). 

In the other limit t±/tu < 1 of two weakly coupled 
chains the magnetic correlation length along the chains di- 
rection becomes larger. Although no magnetic long range 
order exists, a description of the single particle properties 
in terms of a Hartree Fock spin-density-wave (SDW) pic- 
ture turns out to be reasonably accurate |22| . For both 
q± = and q±_ = 7r a dispersive structure is observed with 
a (SDW-likc) gap ~ U separating the photoemission and 
inverse photoemission energy regions. It is worth noticing 
that the low energy electron or hole excited states now 
occurs at momentum q = tt/2 in contrast to the large 
t±/tu limit where they occur at momenta q = 7r (lu < fi) 
and q = (w > fi). The Hartree-Fock treatment correctly 
predicts a bandwidth of order Jii due to the magnetic scat- 
tering (similar to the spinon-like excitations of the single 
chain fl23j). However, it fails to reproduce the broad inco- 
herent background reminiscent of the holon excitations of 
the single chain p3| . 

Away from half- filling (n=l) QMC simulations face 
the well known "minus sign" problem (especially at low 
temperature and large U) which increases the statistical 
errors and, hence, reduces the accuracy of the analytic 
continuation to the real frequency axis. Thus far, QMC 
studies of the doped Hubbard ladder have been restricted 
to U < 4 in the range 1.4 < t±/tn < 2 and for tempera- 
tures larger thani||/8 ^ij. Density matrix renormalization 
group techniques, on the other hand, can currently only 



provide information about static correlations p4| . A re- 
cently developed variational technique based on the use 
of a reduced Hilbert space once the ladder problem is ex- 
pressed in the rung-basis [^5| can produce accurate dy- 
namical results on 2 x 20 clusters at zero temperature and 
finite hole density prj| ]. However, this technique has been 
applied only to isotropric ladders thus far. In the present 
work, alternative approaches have been used. First, follow- 
ing Ref. |22|, a simple estimation of A(q,uj) in the single 
rung approximation has been carried out. This calculation 
is valid in the limit of small i|| and it is useful in order to 
discuss the possible existence of metal-insulator transition 
at commensurate densities such as n = 0.5 or n = 0.75. 
This simple analytical scheme gives also a basis for un- 
derstanding more elaborate numerical calculations. In a 
second step, exact diagonalization studies based on the 
Lanczos algorithm were performed to investigate a broad 
region of parameter space. 



2.2 Local rung approximation: metal-insulator 
transition 

Let us consider the limit where i|| is the smallest energy 
scale, i.e. t\\ <C t± and t\\ <C U. First, A(q±,u) can be 
calculated straightforwardly at densities n — 1 and n = 
0.5 by diagonalizing exactly the single rung Hamiltonian 
for 0, 1, 2 and 3 particles (see Ref. [^2| for details). At 
half-filling one obtains, 

A(0, oj) = a 2 5{uj - f2(2, 1)) + (1 - a 2 ) S(cu - 12(3', 2)) , 
A(ir, u) = (l- a 2 ) S(w - Q{2, l')) + a 2 S(uj - 12(3, 2)X6) 



where a 2 = ^(1 - 



= ) and /2(n, m) correspond to 



^C 2 + (4t_L)2 

the excitation energies of the various allowed transitions 
between a state with m particles to a state with n parti- 
cles. Here, n, n' , n", etc... index the GS and the excited 
states with n particles on a rung of increasing energy. The 
poles of the spectral functions are given, also for increasing 
energies, by 



12(2,1') = 




Q(2,l) = 


> 


f2(3,2) = 


> 


Q(3',2) = 


> 



(U - + (4^)2) + t± , 



(7) 



(U + y/U* + (4*IF) - t± , 



(U + y/U 2 + (4t ± )2) + t± 



The chemical potential /i lies between £2(2, 1) and i?(3, 2) 
leading to the same integrated spectral weight (= 1) in 
the photoemission and inverse photoemission parts of the 
spectrum for all values of U . Hence, as expected, the sys- 
tem is an insulator with a single particle gap of A e h = 

Vu 2 



16t 2 ± 



2t±. However, since the weight a 2 varies 
strongly with the ratio U / At± , the distribution of spectral 
weight changes qualitatively for increasing U as shown in 
Figs. @(a),(b). At small U, a 2 - 1 - j(j^) 2 and one re- 
covers, as in the non-interacting limit, two highly weighted 
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bonding (at an energy around —t±) and antibonding (at 
an energy around t±) bands. When U — > oo, a 2 ^ 1/2 
and, thus, with increasing U spectral weight is transferred 
to bonding and antibonding states further away from the 
chemical potential. In the large U/4t± regime, the sys- 
tem becomes a four-band insulator with 4 (almost) equally 
weighted poles and a Hubbard gap of order U separating 
bonding or anti-bonding states at u> < \x and u> > fx. 

Note that a similar transition from a two-band to a 
four-band insulator has also been observed in QMC stud- 
ies of the half-filled Hubbard ladder at finite f|| and 
fixed value U/t\\ = 8 by decreasing the ratio *j_A|| from 
2 to 0.5. In fact, U/tu = 8 and t±/tu = 2 correspond to 
the intermediate regime U/4t±_ = 1 where, according to 
the previous tii — * estimate, only ~ 15% of the spec- 
tral weight is located in the side bands. For smaller t±, 
more weight appears at the position of these two addi- 
tional structures leading to four bands. In general, for 
arbitrary ratio tj_/t\\, one expects a transition to a four- 
band insulator when U becomes sufficiently large com- 
pared to the largest of the two hopping matrix elements 
i.e. U > Max{£_L,£||}. 

Let us now turn to the discussion of the quarter-filled 
case n = 0.5 where a similar local rung calculation leads 
to, 



A(0, u) = l - 5(u - 0)) + l -a 2 8{lo - Q{2, 1)) 

+ i(l-a a )d( W -/2(2'",l)), 
A(tt, u) = 2 8(uj - f2(2', 1)) +ls(u- 0(2", 1)) , 
where the new energy poles are given by 

fl(i,o) = -tL, 

fl(2 , ,l)=t x , 
f2(2", 1) =t x + U, 
n{2'",l) = Q(3',2). 



(8) 



(9) 
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Fig. 1. Schematic representation of the single particle spec- 
tral function vs frequency, in the t« — > limit. The posi- 
tion of the single particle energy poles are indicated by ar- 
rows whose lengths are proportional to the spectral weights 
associated to the corresponding transitions. Arrows pointing 
upwards (downwards) correspond to q± — (q± = n). The 
photoemission peaks (occupied states) and the inverse pho- 
toemission peaks (empty states) correspond to full line and 
dashed line arrows, respectively. The spectra are shown for ra- 
tios ^ ~ 1/4 ((a), (b) and (c)) and £ ~ 2 ((d), (e) and 
(f)) and for electron densities n = 1, n = 0.5 and n = 0.75 as 
indicated. 



Since the chemical potential is located exactly between 
J7(1,Q) and 17(2', 1), the system is an insulator for all 
values of U and (sufficiently) small tu (compared to U). 
However, A(g±,ui) exhibits completely different forms at 
small and large U couplings as observed in Figs. |l](c),(d). 
At small U, as in the non-interacting case, the bonding 
states and antibonding states are separated by an en- 
ergy of order 2t±. However, each structure is split by a 
small energy proportional to U and the chemical poten- 
tial lies between the two q± — sub-bands. For large U, 
the gap becomes of order 2t± and upper Hubbard bands 
(of almost equal weights) of the bonding and anti-bonding 
states appear at an energy ~ U higher. Although this pic- 
ture does not take into account tn, the role of a small charge mode) is stable for arbitrary small tn. At n = 0.75 



discussed in the next section. Since the single particle ex- 
citation gap A eh = \(U - ^U 2 + 16t 2 ± ) + 2t± is of the 
order of the smallest of the two energy scales U/2 and 
2t±, the insulating phase is then restricted to the range 
4i|| < Min{t7/2,2*x}. 

The existence of a metal-insulator transition is, in fact, 
specific to quarter filling (besides the half-filled case which 
is always insulating). A simple argument is here presented 
to show that at other (commensurate) densities such as 
n = 0.75 the metallic phase (i.e. with at least one gapless 



*H can be easily discussed qualitatively. In fact f|| is ex- 
pected to give a dispersion in the chain direction and a 
width to the various structures discussed here. When 4*ii 
becomes comparable to the single particle excitation gap 
A e h , bands will start to overlap and a transition from the 
insulator to a metallic state is expected fl27||, as will be 



a local rung calculation of A(qj_,u) requires to consider 
as a GS two decoupled rungs on 4 sites with 2 and 1 
particle, respectively. The spectral function is then given 
straightforwardly by the average of the spectral function 
Eqs.(||) and (§) at densities n = 1 and n — 0.5. However, 
the location of the chemical potential is a subtle issue: 
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since the states at the energy lo — 0(2, 1) are completely 
filled (empty) for n = 1 (n = 0.5), it is clear that this 
state will become partially filled at n = 0.75 so that the 
chemical potential is pinned at this energy. Consideration 
of the spectral weights of the excitations shows immedi- 
ately that, for arbitrary small tn, the band centered at 
lo = [2(2, 1) (of weight |a 2 ) is always |-filled leading to 
a metallic behavior. In this case, an additional interac- 
tion, e.g. between nearest neighbor sites along the chains, 
would be required to produce a metal-insulator transition. 
Schematic representations of A(q±,Lo) at n = 0.75 are 
shown in Figs. [j](c,f) at small and large U. At small [/, 
as expected, the two bonding and anti-bonding structures 
separated by ~ 2t± are clearly visible and the bonding 
states at the lower energies are partially occupied. In this 
limit, U leads essentially to small splittings of the various 
structures into sub-bands (as for n = 0.5). For large U, 
the spectral function is qualitatively very different with 2 
distinct bands around — tj_ and t± for both q± = and 
q± = 7r states. However, the upper Hubbard band around 
an energy ~ U is formed of two peaks (separated by 2t±) 
for q± = 7r while only one peak is present for q±_ = 0. 

Finally in this section a brief discussion of the case of 
the t — J ladder is included. Since this model describes 
only the low energy properties of the Hubbard model, the 
corresponding spectral functions in the iii — > limit can 
be obtained easily from the previous ones by discarding 
the high energy peaks whose energy scales as U for large- 
[/, setting a 2 = 1/2 and expanding energies to first order 
in J± = At\/U . In fact, it can be easily shown that the 
same expressions hold for the t — J model (with arbitrary 
J±). Note that, due to the projection of the high energy 
states, the spectral function of the t— J model follows the 
new sum-rule J A(q±,Lo) dtu — it 2 (instead of 1), where 
x = 1 — n is the doping fraction. At half-filling one gets: 

4(7r,w) = i*(w-(-t x -Jx)), (10) 

with the chemical potential located at an higher energy 
(~ U/2). Similarly, at quarter-filling one obtains, 

A(0, w) = \ 8(oo - (-i x )) + \s(u - (t x - J X )) , 
A(tt,w) = ^(w-*x), (11) 

with the chemical potential located between —t± and tj_ — 
J±. It is interesting to notice that, when J± exceeds 2t±, 
the electron-like excitation becomes lower in energy than 
the hole-like excitation. This signals the onset of phase 
separation or, alternatively, some sort of charge localiza- 
tion/ordering (such as charge density wave ordering). Phys- 
ically, this occurs when the magnetic energy gain of a sin- 
glet on a single rung becomes larger than the kinetic en- 
ergy of two particles on individual rungs. 



2.3 Exact diagonalization results: Hubbard model 

Let us now investigate the dynamical properties of the 
Hubbard and t — J models for arbitrary parameters us- 
ing exact diagonalization techniques. Cyclic 2 x L ladders 
are diagonalized and the (zero temperature) particle spec- 
tral function is obtained exactly by a standard continued- 
fraction procedure. Although in practice one is limited to 
L = 8 (for the Hubbard model), both periodic (PBC) and 
anti-periodic (ABC) boundary conditions can be used to 
consider a sufficiently large number of momenta q\\ = nj-, 
rc=0,2J-l. 

The case of the Hubbard ladder will be considered first, 
before focusing on the low energy excitations described by 
the t — J model. The spectral function A(q, lo) at a density 
of n — 0.75 is shown in Figs. |](a,b,c) for U = 8 and several 
values of t± ranging from 2.5 down to 0.5. Note that both 
PBC and ABC have been used in Fig. ||(b) while, in order 
to reduce CPU time, only ABC (PBC) have been used 
in Fig. ||(a) (Fig. ||(c)). Two sharp structures separated 
by an energy proportional to t± can be attributed to a 
bonding and an anti-bonding band. At the largest ratio 
of t±/tu = 2.5 that have been considered, the spectrum 
exhibits some features of Fig. [l](c) obtained in the local 
rung approximation at small coupling: (i) in the photoe- 
mission part, a q± — sub-band of small spectral weight 
can be observed at an energy of about U/2 from the main 
q± = band crossing the chemical potential; (ii) & q± = n 
upper Hubbard band appears at an energy oc U away from 
the main (empty) q± = n band. On the other hand, some 
tiny structures characteristic of the strong coupling limit 
(Fig. 0(f)) can also be observed: (i) a small spectral weight 
exists at lo < /i (around to ~ —5) for q± = ir together with 
(ii) a quite small q± = upper Hubbard band at lo > /i. 
Interestingly enough, these features become more impor- 
tant for t± = 1.5 as shown in Fig. ^(b) which corresponds, 
in fact, to a larger ratio U /At± ~ 1.3. 

With decreasing electron density, the respective posi- 
tion of the two main bands and the position of the Fermi 
level seems to evolve as in a rigid-band scheme. How- 
ever, there are important differences: (i) the bandwidth is 
strongly reduced specially at smaller t±/tu ; (ii) the excita- 
tions become sharper when the band crosses the chemical 
potential. To the best of our knowledge, this is the first 
observation in a numerical study of the broadening of the 
"quasi-particle" -like peaks excitations as one moves away 
from the chemical potential. 

For a larger hole doping and working at a commensu- 
rate value of n = 0.5 qualitative changes can take place in 
the spectral function at sufficiently large U and t±. Data 
are shown in Fig. |3|. For t± — 0.5 the two partially filled 
bonding and antibonding bands can be observed together 
with their corresponding upper Hubbard bands at higher 
energy. As expected from the previous fii — > analysis, 
the spectral weight of the q± = upper Hubbard band, at 
fixed U, gets strongly reduced for increasing i.e. for a 
decreasing ratio U/t±. At large enough t± , a gap appears 
in the q± = structure, leading to two sub-bands and an 
insulating behavior in agreement with the local rung cal- 
culation. Such a metal-insulator transition is induced by 
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Fig. 2. Spectral function A(q,io) of the Hubbard ladder for 
[7 = 8 and n = 0.75. The left and right sides correspond to the 
bonding (k y — 0) and anti-bonding states (k y = 7r), respec- 
tively, and fez runs from to tv from the top to the bottom. 
The position of the chemical potential is indicated by a vertical 
dotted line, (a), (b) and (c) correspond to t± = 2.5, t± = 1.5 
and t± — 0.5, respectively. 



Fig. 3. Spectral function A(q,tu) of the Hubbard ladder for 
U = 10 at quarter-filling n — 0.5. The left and right sides 
correspond to the bonding (k y — 0) and anti-bonding states 
(k y = 7r), respectively, and k x runs from to tt from the top to 
the bottom, (a), (b) and (c) correspond to U = 10 with t± = 5, 
t± = 2.5 and t± = 0.5 respectively. 
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a combined effect of t± and U : when t± is large enough 
the lower band becomes half-filled and a finite U, leading 
to relevant Umklapp scattering, can then open a gap. 

2.4 Exact diagonalization results: t — J model 

In order to study more precisely the influence of doping at 
small energy scales let us now focus on the t— J ladder |]28| . 
Results at small hole densities n — 0.875 and n = 0.75 are 
shown in Figs. |^ and Figs, [s] and are consistent with the 
previous results on the Hubbard model. 

Let us first discuss the role of the hole doping x = 
1 — n, for the largest value of t± = 2 considered here 
(see Fig. [|(a) and Fig. ||(a)). For this choice of param- 
eters, the x-dependence can be qualitatively understood 
from the single rung picture. For U — > the GS contains 
a density of 2x singly occupied bonds and 1 — 2x dou- 
bly occupied bonds. By combining the spectral functions 
at n = 0.5 and n — 1 with the respective weights, one 
obtains a simple picture of the influence of doping con- 
sistent with the numerical results at small (but finite) tn. 
The q± = main structure (which is the closest to the 
chemical potential at half-filling) becomes partially filled 
with a weight of x/2 in the inverse photoemission part 
oj > fx. Note that the dispersion of the band is especially 
flat in the vicinity of the chemical potential at small x. 
With increasing doping, weight is transferred from this 
structure (of total weight 1/2 — x/2) and from the up- 
per Hubbard band (not described by the t — J model) to 
q± = states further away from the chemical potential. 
This leads to an emerging structure of weight x at an en- 
ergy of ~ 2t± — J± below the main band. Physically, in a 
photoemission experiment, these small peaks correspond 
to processes where an electron on a singly occupied rung 
is removed by a photon and leaves behind an empty rung. 
Note that this structure becomes particularly strong at 
quarter filling (as seen in Fig. ^(a)) where it carries 1/2 of 
the total spectral weight (normalized to 1). In the q±_ = n 
sector, the main structure in the photoemission part of the 
spectrum u> < fi (barely seen in the case of the Hubbard 
model for the parameters chosen in the previous study) 
is also loosing spectral weight upon doping with a total 
weight of 1/2 — x. The missing weight (and some addi- 
tional spectral weight from the upper Hubbard band) is 
transferred into the inverse photoemission spectrum lead- 
ing to an emerging band of total weight 3x/2 at lj > /i. 
Such states, obtained by suddenly adding an electron on 
a singly occupied rung could be seen in an inverse photoe- 
mission experiment. At quarter filling n = 0.5, as seen in 
Fig. [|(a), the transfer of spectral weight is complete and 
the u> < fi, q±_ = it structure has totally disappeared. 

At smaller values of t± (see Figs. ||(b,c) and Figs. ||(b,c)) 
the two separate structures, both for q± = or q± = it, 
merge into a single broad structure. The data can be fairly 
well described by (i) q± = and q± = it bands dispersing 
through the chemical potential and (ii) a broad incoher- 
ent background extending further away from the chemical 
potential towards negative energies. Note that, similarly 
to the previous case of the Hubbard model, the peaks of 
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the band-like feature seem to become narrower when they 
cross the chemical potential as expected in a Fermi liquid 
description. 

This Section ends with a short discussion on the pos- 
sible existence of small single particle gaps in the previous 
data. The quarter-filled case is qualitatively different from 
the low doping regime. At n = 0.5 the tii — > analysis un- 
ambiguously predicts the existence of a gap in the single 
particle spectrum at sufficiently large t± and U. However, 
when the q± = it structure is not completely empty (i.e. 
totally located in the u> > fi region of the spectrum), as it 
is the case in Figs. ||(c) and ||(b,c), no gap is expected as 
it is clear in the numerical data. Then, a metal-insulator 
transition is expected by increasing t± but whether this 
transition is driven by t±_ alone is still unclear. The data 
of Figs. §(a,b) corresponding to a situation where the an- 
tibonding band is clearly unoccupied do not allow to ac- 
curately determine a critical value of t± at which the gap 
starts to grow. However, we have checked numerically (not 
shown) that, by reducing t\\, the spectrum of Fig. |3](a) 
smoothly evolves into the spectrum obtained above in the 
single rung approximation, e.g. exhibiting a well defined 
gap at the chemical potential. 

At small doping, on the other hand, the physical origin 
of a small single particle gap would be quite different. In 
this case, it would be related to the formation of pairs. In 
the t|| — > limit, pairs become stable only when J± > 2t±, 
i.e. when the magnetic energy on a rung exceeds the ki- 
netic energy loss. Otherwise, for J± < 2t±, the spin gap 
is immediately destroyed by doping (strictly for tn = 0) 
since the presence of singly occupied rungs leads to new 
low-energy spin-1 excitations in the n = 1 spin gap (of 
order Jj_). Therefore, intermediate ratios of tj_/tu seem to 
be more favorable for pair binding. Although spectra like 
those shown in Figs. [|(b) and Figs. ||(b) are not incon- 
sistent with the presence of a small gap at the chemical 
potential, the study of pair binding from an investigation 
of the spectral function at small energy scales around the 
chemical potential is a difficult task. In order to clarify this 
issue, a complementary study of static physical quantities 
is shown in the next Section. 



3 Superconducting properties 

3.1 Pair binding energy 

In the limit where J± is the largest energy scale, forma- 
tion of hole pairs are favored on the rungs in order to min- 
imize the magnetic energy cost. In fact, this simple naive 
argument breaks down when tj_ > Ji/2 since holes on 
separate rungs can then benefit from a derealization on 
each rung. In the large t± limit, a simple 4-sites (2 rungs) 
calculation shows that for J±/2 << t±, the pair binding 
energy (which, as defined below, should be positive if a 
bound state exists) behaves as J?JAt±_ — 2t\\. Very tightly 
bound hole pairs are then not stable in the intermediate 
regime. However, extended pairs have been shown to exist 
in some regime [[5|]6lfj9|| . In this Section, the stability of 
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respectively. 



the hole pairs as a function of the anisotropy ratio will be 
investigated in detail. 

To study the onset of pair binding, irrespective of the 
actual size of the pair, the two- hole binding energy defined 

by 

A B {L) = 2E(L, 1) - E(L, 2) - E(L, 0) , (12) 

is considered, with E(L, Nh) being the GS energy of the 
2 x L ladder doped with N h holes. E(L, N h ) (N h = 0, 1, 2) 
have been calculated out— J ladders with periodic or anti- 
periodic boundary conditions along the legs direction and 
considering sizes up to L = 13. As previously, it has been 
assumed for convenience that J± = J||(tx/t||) 2 - Typical 
finite size behaviors of Ag(L) are shown in Figs. 0(a,b) 
for different anisotropy ratios. Some caution is obviously 
needed in order to extrapolate the results to the ther- 
modynamic limit. At large t± for a fixed choice of the 
boundary conditions very regular oscillations in As vs L 
appear (as observed for instance in Fig. 0(a)). However, 
the data set corresponding to ladders with an even (odd) 
number of rungs (L) and PBC can be combined with the 
data set corresponding to ladders with an odd (even) num- 
ber of rungs and ABC. The two resulting curves exhibit a 
smooth behavior (full lines) which allow an accurate ex- 
trapolation to the thermodynamic limit. For intermediate 
values of t±/t ~ 1.25, the two data sets merge into a sin- 
gle curve and finite size corrections become particularly 
small. 

For smaller ratios such as £j_/t|| < 1 (see Fig. 0(b)), 
there is a qualitative change of behavior. In this case, 
the data corresponding to ladders with an odd number of 
rungs (L = 2p + 1) have to be distinguished from the data 
obtained for L = 2p. Indeed, in this parameter regime, 
the spin correlation length along the chains becomes of 
the order of the system size so that ladders with L = 2p 
and 2 holes suffer from a small magnetic frustration in- 
duced by the boundary conditions. An accurate finite size 
scaling analysis can nevertheless be realized by consider- 
ing the data for L = 2p + 1 which show again a very 
systematic behavior as a function of L; in fact, although 
the behavior of the data sets corresponding to ladders with 
L = 4p+l (L — 4p+3) rungs and PBC follow only roughly 
the same trend as the data set corresponding to ladders 
with L = Ap + 3 (L — Ap + 1) rungs and ABC, averaging 
over PBC and ABC leads to a single remarkably smooth 
behavior of Ab vs L as observed in Fig. 0(b). A simi- 
lar procedure can be followed for the data obtained with 
L = 2p as shown also in Fig. 0(b) but this extrapolation is 
probably less reliable for the reasons stated above. In all 
cases, the L — > oo extrapolation is performed according to 
A B {L) = A%° + A\ exp (-L/0, where A™, £ and A are 
free parameters determined from a fit to the data. 

The extrapolated values of Ab are displayed in Fig. 
as a function of t±_ for Jii = 0.5. A positive binding en- 
ergy (implying the stability of the hole pair) is obtained 
for all parameters considered here. Recent DMRG work 
using clusters with up to 2 x 30 sites and two holes have 
calculated the binding energy in the isotropic limit [ po| . 
The result Ab ~ 0.2 reported there is very close to the 
number found by our size extrapolation in Fig. 0, giving 
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extra support to our procedure. For comparison, the spin 
gaps Aq and A2 in the undoped GS and 2 hole doped GS, 
respectively are also shown. The scaling behavior of A% is 
very similar to that of Ab as can be seen in Figs. @(c,d) 
(although finite size effects are larger) and the same pro- 
cedures to obtain the extrapolations to L = oo have been 
used. Our extrapolations for Aq are in good agreement 
with previous ED J3l| and QMC |32| estimates or with 
the weak coupling limit behavior Aq ~ 0A1J±. The be- 
haviors of As and A2 with t± are very similar, and both 
have a pronounced maximum of Ab at tx/t ~ 1.25. In 
fact it is expected that in the 2 hole doped GS the lowest 
S = 1 excitations can be obtained by either (i) break- 
ing up a hole pair and flipping one of the unpaired spins 
or (ii) making a spin excitation away from the hole pair 
(which is supposed to have a finite size). Therefore, one 
expects A 2 — Min {Ab, A }. Our data indeed suggest two 
regimes: (i) for tx/tn < 1.25, the binding energy of a pair 
is larger than the spin gap in the undoped system. This 
extra stability of the hole pair is probably due to the 
strong antifcrromagnetic correlations within each chain. 
Note that, in this regime, there are some small discrepan- 
cies between A2 and Aq. This effect is probably due to the 
fact that the size of the pair is particularly large for such 
parameters and the two regions corresponding to the un- 
doped spin liquid and the hole pair cannot truly separate 
in the clusters that can be handled computationally lead- 
ing to strong finite size effects for A 2 - (ii) for tx/t\\ > 1.25, 
the hole binding energy strongly decreases and the lowest 
spin excitation is obtained by breaking a hole pair i.e. 
A2 ~ Ab- Note that for even larger ratios tx/tn (typ- 
ically ix/iii > 2.5), the binding energy increases again. 
Clearly, this is an artificial effect due to the fact that, 
in our model, the rung magnetic coupling scales like t\ 
and becomes unphysically large compared to tj_ for large 
enough t±. In that case, Ab — J± — 2t± which approaches 
the spin gap Aq for large t±. 

It is interesting to compare the results of Fig. || with 
the previous study of the collective modes of the t — J lad- 
der flipf . On general grounds, two collective spin modes of 
momenta q± — and q± — tt are expected in a doped spin 
ladder. Both modes are gapped at moderate doping ]l0| . 
From a careful examination of the quantum numbers of 
the various spin excitations shown in Fig. |[ one can safely 
study, at vanishing doping (i.e. for 2 holes in an infinitely 
large system), each low energy excitation. The collective 
qx = 7r spin mode corresponds to the spin excitation of 
energy Aq characteristic of the undoped system (crudely 
an excitation of a singlet rung into a triplet) . On the other 
hand, the qx = spin mode is associated to the breaking 
of a hole pair of energy Ab- From our previous analy- 
sis of the data, a level crossing occurs between these two 
types of excitations around t±_ ~ 1.25 producing a cusp- 
like maximum of A^- 

Materials corresponding to the regime tx/t\\ > 1.25 
should be particularly interesting to be studied by Inelas- 
tic Neutron Scattering (INS) experiments at small doping. 
Indeed, the above calculation predicts that, under light 
doping, spectral weight in the dynamical spin structure 




Fig. 7. Finite size scaling behaviors as a function of the inverse 
of the ladder length for Jy = 0.5. Filled circles (open squares) 
correspond to PBC (ABC). The values of t± are shown on the 
plot. The full lines correspond to the finite size scaling laws 
used for the extrapolations to L = 00. (a) Two hole binding 
energy for t± = 2.25; (b) Two hole binding energy for t± = 
I/a/2; (c) Finite size behavior of the triplet gap in the GS with 
2 holes for t± — 2.25; (d) Finite size behavior of the triplet 
gap in the GS with 2 holes for t± — 1/V2; In (b) and (d), 
the triangles correspond to averages between the PBC and the 
ABC data and the sizes of the error bars correspond to the 
absolute value of the difference between the two sets. Open 
(filled) symbols correspond here to L odd (even). 
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Fig. 8. Extrapolated two hole binding energy (open symbols) 
as a function of the anisotropy r a = t±/tn for Jii = 0.5. Spin 
gaps in the half-filled and two hole doped GS are also shown 
for comparison (filled symbols). 
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factor ^(q, ui) should appear within the spin gap of the 
undoped material. This new q± = magnetic structure 
whose total weight should roughly scale with the doping 
fraction corresponds to the excitations of hole pairs into 
two separate holes in a triplet state. The corresponding 
energy scale for such excitations can be much lower than 
the spin gap of the undoped spin liquid GS (q± = n) . 

The maximum observed in As for the t — J model at 
Jii = 0.5 as a function of t± has similarities with the be- 
havior of the pair-pair correlation obtained in the Hubbard 
model at very small hole doping |24| n = 0.9375 which 
also shows a maximum (around t± ~ 1.4 for U = 8). In 
Ref. pi! , this particular value of t± was associated with 
the situation where the chemical potential coincides al- 
most exactly with the top of the lower bonding band and 
with the bottom of the upper antibonding band. In that 
case, one expect a particularly large density of state at 
the chemical potential (see also Ref. j36|). However, such a 
correspondence was made possible at smaller U only (due 
to difficulties to obtain accurate QMC calculations of dy- 
namical quantities at intermediate and large values of U). 
The spectral function A(q, ui) shown in Fig. [I|(b) was ob- 
tained in the two hole GS of the 2x8 ladder for a choice 
of parameters (Jii = 0.4 and t± = 1) close to the ones 
producing the maximum of A B in Fig. |. Fig. |(b) clearly 
shows a large density of states in the vicinity of the chem- 
ical potential due to the flatness of the dispersion around 
q = (0, 7r) or q = (ir, 0). This situation corresponds to the 
cross-over between the two band and four band insulator 
regimes observed at half- filling It is also interesting 
to note that a small depression of the density of state is 
visible in Fig. ||(b) at the chemical potential. This could 
be interpreted as a small gap associated to the existence 
of a bound pair. More generally, in the so called C1S0 
phase [ ^|Jll| ,16|] where the spin gap survives, one expects 
to see its signature in A(q, ui) as a gap at the chemical po- 
tential. However, the energy scale of the spin gap is small 
(see e.g. the order of magnitude of A^ in Fig. |J) com- 
pared to the various features that appear in A(q, ui) and 
thus, in most cases, its manifestation in A(q, ui) cannot be 
observed on small lattices. In the recent studies using a re- 
duced Hilbert space, the observation of a gap caused by 
pairing in the spectral function required the use of clusters 
with 2 x 16 and 2 x 20 sites 



3.2 Pair-pair correlations 

ED studies supplemented by conformal invariance argu- 
ments suggest that in the doped spin gap phase (C1S0) of 
the isotropic t — J ladder (where pairs are formed accord- 
ing to the previous analysis) algebraic superconducting 
and 4/cf-CDW correlations are competing Q. At small 
J/t ratio, the CDW correlations dominate while above 
a moderate critical value of J/t coherent hopping of the 
pairs takes over. The aim of the present Subsection is to 
investigate the role of the anisotropy t±/t\\ by a direct 
calculation of the pair-pair correlation as a function of 
distance. As previously, in the case of the t — J model, a 
rung magnetic coupling J± = J\\(tj_/t\\) 2 is used. 



Superconducting correlations can be evidenced from a 
study of the long distance behavior of the pair hopping 
correlation, 



C s {r-r') = (A\r)A{r')), 



(13) 



where A'{r) is a creation operator of a pair centered at 
position labeled by r. Although the best choice of A\r) 
clearly depends on the internal structure of the hole pair [B3| 
as discussed later, it should exhibit general symmetry prop- 
erties associated to the quantum numbers of the hole pair 
found in the previous Subsection: (i) A^ (r) is a singlet op- 
erator and (ii) it is even with respect to the two reflection 
symmetries along and perpendicular to the ladder direc- 
tion (and centered at position r). The static correlation 
function of Eq. ^| can be interpreted as a coherent hop- 
ping of a pair centered at position r to a new position 
r'. 

According to conformal invariance, in a strictly ID lad- 
der (which is the case studied here) the pair hopping cor- 
relation exhibits a power-law behavior at large distances 



C S (r 



1/r- 



(14) 



where the exponent K p was calculated in the weak cou- 
pling limit |lT| or in the isotropic t — J ladder by ED 
methods using conformal invariance relations ||l3f . Super- 
conducting correlations dominate when K p > 1/2 which 
occurs for J/t > 0.3 in the lightly doped isotropic t — J 
ladder |l3f. Using a DMRG approach, the behavior of 
Cs{f — T) with the usual BCS bond pair operator, 



(15) 



can also be obtained directly, leading, in the case of the 
isotropic t — J model (l2j , to a good agreement with the 
ED results. More recently, this study was extended to the 
anisotropic Hubbard ladder Q showing a pronounced 
peak of the long-distance pair correlations as a function 
o£t±. 

Here, as a complementary study of the analysis pre- 
sented for the binding energy in the previous Subsection, 
the behavior of the pair correlation function of the BCS- 
like operator of Eq.( [l5|) is compared against the case 
where a spatially-extended pair operator is used. The first 
motivation to introduce this new pair operator is due to 
the structure of the hole pair; indeed, it turns out that 
configurations in which the two holes sit along the diago- 
nal of a plaquette carry a particularly large weight in the 
2-hole GS both in the case of the 2D t — J model or in 
the case of the t — J ladder |29j . This feature seems coun- 
terintuitive in a two-hole bound state of d x 2_ y 2 character, 
as it is the case e.g. in 2D (for ladders, this symmetry is 
only approximate), since the pair state is odd with respect 
to a reflection along the plaquette diagonals. However, it 
has been observed [Q that retardation provides in fact a 
simple physical explanation of this apparent paradox. Sec- 
ondly, it is clear that pairs extending into a larger region 
of space can acquire more internal kinetic energy and they 
are less sensitive to short distance electrostatic repulsion. 
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Fig. 9. Pair-pair correlation function vs distance calculated 
on 2 x 8 clusters at density n = 0.75 with PBC in the chain 
direction. The values of the anisotropy r a = tx/t\\ are indicated 
on the plot, (a) rung-rung correlations in the Hubbard ladder 
for U = 10; (b) rung-rung correlations in the t — J ladder for 
J|l = 0.4; (c) plaquette-plaquette correlations (open symbols) 
in the t — J ladder for Jn — 0.4. For comparison, some of the 
correlations of the rung pair operator of (b) are also reproduced 
(small full symbols) on the same plot. 



To study the influence of the spatial extension of the 
pair operator A{r) on the pair-pair correlation, following 
Ref. [p3| here a plaquette pair operator is defined as 

A(i + 1/2) = (Si,2 - S<+l,l) • T,.|:,. |.2 

— — Sj+i^) • Tj+iji;,^ (16) 

where Ti, a - 3 ,(3 = \ci, a - a (a y cr) aa >Cj^- ta > is the regular (ori- 
ented) spin triplet pair operator |35|] , Physically, A^(i + 
1/2) creates a singlet pair centered on a plaquette in a 
d x 2_ y 2 state with holes located along the diagonals of the 

plaquette (at distance y/2) . The interpretation of this op- 
erator is simple: starting from a hole pair located on a 
rung, the hopping of one of the holes by one site along 
the leg-ladder leaves behind a spin with the opposite ori- 
entation than the local AF pattern. This argument nat- 
urally leads to a 3-body problem Q involving a triplet 
hole pair and a local spin flip (of triplet character). For- 
mally, this picture is equivalent to introducing some re- 
tardation in the usual BCS operator i.e. the two holes can 
be created at two different times separated by an amount 
r e.g. by applying c^ijf (r/2)c ij2; x (— r/2) on the AF back- 
ground. The expansion of this new operator to order r 2 
then leads to the various terms of Eq. ( [l6]) . Alternatively, 
A<(i + 1/2) can also be viewed as the simplest d x 2_ y 2 op- 
erator of global singlet character creating a pair on the 
diagonals of a plaquette. This result can be deduced from 
simple symmetry considerations p3| . 

Our results for C s (r) in the case of the rung BCS- like 
operator are shown in Fig. ^(a) and (b) for the Hubbard 
and t — J ladders, respectively. Both sets of data are con- 
sistent with the power law decay and show a clear increase 
of the correlations at intermediate distances. In the case 
of the t — J ladder at n = 0.75, the maximum occurs for 
t± ~ 1.5, a value slightly larger than the characteristic 
value corresponding to the maximum of Ab- According 
to Figs. H and [5] showing the single particle spectral func- 
tions for almost identical parameters, this specific value 
of t± seems to correspond to the case where the chemical 
potential sits in the vicinity of a maximum of the density 
of states generated by very flat bands at the band edge (as 
suggested in Ref. |24j and in agreement with the general 
ideas discussed in Ref. |56|). On the other hand, it is likely 
that the maximum of As does not occur at exactly the 
same value of t± but rather at a somewhat smaller value. 

The plaquette pair-pair correlations are shown in Fig. ^(c). 
At short distance r = 1, the correlations are suppressed 
reflecting the spatial extension of the pair operator. At 
larger distances, r > 2, a significant overall increase is ob- 
served compared to the case of the rung operator, showing 
that indeed the use of "extended" operators to capture the 
usually weak signals of superconductivity in doped anti- 
ferromagnetic systems is a promising strategy |3j. Note 
that, apart from this overall factor, the functional form of 
the decay seems to be identical to the one obtained for 
the rung operator (as can be checked quantitatively). 
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4 Conclusions 

In this paper dynamical properties of anisotropic ladders 
have been investigated using the one-band Hubbard and 
t — J models. An analysis based on the local-rung ap- 
proximation explains a considerable part of the numerical 
results. In particular, the existence of a metal-insulator 
transition at quarter filling which can be justified in such 
an analysis was indeed numerically seen for increasing 
anisotropy ratio. Flat quasiparticle dispersions at the chem- 
ical potential are observed in regions of parameter space 
where pairing correlations are robust. A finite-size scaling 
of the binding energy and the spin-gap show that these 
quantities change with the anisotropy ratio in a manner 
similar as the pair correlations do. In agreement with pre- 
vious results, it is observed that superconducting corre- 
lations are maximized for anisotropic systems, with cou- 
plings along rungs slightly larger than along the legs. 
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